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FOREWORD 


Thin  research  was  support ed  by  the  Advanced  Research 
Projects  Agency  of  the  Department  of'  Defense  and  v/as  moni¬ 
tored  by  the  Air  Force  Office  of  Scientific  Research  under 
Grant  AP0SR-72-2305. 


ABSTRACT 


The  effectiveness  of  large  arrays  in  mapping  the 
noise  field  is  well  known .  This  paper  describes  an  at¬ 
tempt  to  estimate  the  noise  fields  by  using  a  small  number 
of  censors.  The  noise  fields  arc  defined  by  their  power 
densities  in  the  frequency-wavenumber  cnacc  and  their 
validity  will  be  judged  by  comparing  coherence  estimates 
of  real  data  with  coin a  .  computations  on  the  basis  of 

the  models.  The  real  data-base  have  been  recordings 
from  Oyer  array  -  the  first  large  installation  in  the 
NOROAR  area. 
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INTRODUCTION 


During  the  last  few  years  methods  fox*  mapping  the 
noise  fields  have  been  presented  in  literature  and  re¬ 
sults  from  different  sites  have  been  presented  in  numer¬ 
ous  reports. 


A  commonly  used  and  excellent;  way  of  presenting  che 
results  is  by  displaying  the  power  density  as  a  function 
of  the  f i eouency  and  the  wavenumber,  thus  giving  the 
distribution  in  asimutli  and  velocity  of  the  noise  fields. 
However,  to  map  the  noise  field  this  way  one  must  have 
access  to  data  from  large  and  properly  spaced  arrays  and 
of  course  the  results  are  valid  only  at  or  near  the 


array  sites. 

Now  the 
large,  and  t 
ally  limited 
one  oi’  a  few 
we  have  .inve 
noise  field 


number  of  large  aperture  arrays  are  not  very 
he  number  of  sample  points  in  space  arc  usu— 


.  The  normal  ease  will  be 
(2-3)  sensors  at  each  site 
stig'ited  the  possibility  of 
in  frequency  --  wavenumber  s 


recordings  from 
.  In  this  paper 
estimating  the 
pace  by  using  the 


experimental  data  i ron  only  a  snail  number  of 
The  procedure  has  boon  to  design  noise  models 
f ic  power  distributions  and  to  check  the  model 
rence  computations. 


sensors . 
with  speci- 
s  by  coho- 
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THE  NOISE  MODEL 


Before  defining  the  noise  model,  let  us  point  out 
a  few  requirements  that  must  he  met:  First  of  all  the 
model  must  be  simple  enough  to  be  rant  hem-  it  i  cully  for¬ 
mulated  and  to  allow  the  calculation  of  the  oararaeter 
of  interest.  Secondly  it  should  not  deviate  too  much 
from  the  noise  fields  as  mapped  by  using  arrays  in  the 
same  area,  and  thirdly  it  must  explain  certain  neculiar 
observations  such  as  variation  of  coherence  with  direct¬ 
ion  (Rygg  and  al.  1969).  With  these  restrictions  in 
mind  we  define  the  model  as  follows:  The  theoretical 
noise  field  consists  of  a  number  of  plane,  uncorrelated 
wavetrains  approaching  from  all  azimuths  and  distributed 
over  a  certain  velocity  range.  (Fig.  1).  Each  wave train 
is  assumed  to  have  a  flat  spectrum  inside  the  frequency 
band  of  interest  for  our  computations  (white  noise),  and 
the  power  density  is  distributed  with  varying  strength 
along  the  periphery. 

The  reason  for  chosina  this  model  instead  of  disc 
noise  sources  or  a  combination  of  disc  noise  and  fixed 
velocity  arc  noise  is  that  we  have  experienced  that  this 
model  in  a  good  approximation  to  the  experimentally  esti¬ 
mated  noise  field  (Bun gum  and  al.  1971).  In  the  node! 
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v.’o  Invc  also  aLlowed  for-  some  now  or  variation  associated 
with  vary  in':  velocity,  thus  taking  into  account  eneryv 
connected  with  different  mode::  of  propagation.  On  evalu- 
catm  •  tlie  thooi-etical  coherence  function  for  thin  model 
we  follow  tile  linen;  of  Murdoch  and  Pfluke  (1970):  fihc 
periphery  i:  divided  into  K  dircrete  directions.  Prom 
each  direction  we  assume  that  J,  discrete  wavetrains  pro¬ 
pagate  with  different  velocities,  carrying  different 
amounts  of  power.  The  total  number  of  ware  trains,  ren- 
chin'-'  a  sensor  is  then  K>  L.  The  output  in  the  transform 
domain  is: 


3(f) 


K.  L 

=  )_  Dm  (f>  •  H,c,i 

k=1  1=1 


Here  Ak>1  is  the  Fourier  transform  at  a  spatial  reference 
point  o4  the  wavct’wiin  with  direction  index  k  and  velocdt 
xncu\  1.  Hj:^(l)  ?.s  o  transfer  function  expressing  the 

cl  feet  of  the  medium  f • om  the  reference  point  to  the  sen¬ 
sor.  (Wo  assume  the  instruments  to  ho  identical). 

The  cross-spectrum  between  two  sensors  (1  and  2) 
is  then: 


>1?(f)  =  s  1  (f )  •  s2(f)  = 


K  L_  K  1. _  _ 

D-  E  E  h,l(f)  •  An,n(f)  • 

k=1  1=1  rn=  1  n-1 


Tho  bars  represent  complex  conjugates .  How,  since  the 
wavetrnins  arc  mutually  uncorrelated : 


E 


•  Vn(f)j 


=  P,7,  ->(f),  when  :n=k  and  n=l 

*'»  ± 

=  0  otherwise. 


Here,  F.V,  (f )  is  tho  (auto) -power  spectrum  of  the  wave- 
train  comin.a  from  the  k'th  direction  and  propa^rt in."  with 
the  velocity  which  in  tied  to  velocity  index  1. 

If  wo  ncylect  attenuation  and  dispersion  across  the 
site,  the  transfer  functions  ±(f )  represent  mer(>ly 
uhnsc  C elo.ys .  Thus,  if  t1  ,  is  the  time  reouired  for 
n  specific  v/avn train  to  nans  from  the  reference  uoint  to 
sensor  1,  the  associated  transfer  function  car.  be  written: 

— r.i)  t  .  i  *1 

Hl,kjl(f)  =  C  l’K’1  (»■«  =  2TTf) 

Then  we  have: 
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II 


i(Jl, i  -i'-’t ,,  ,  , 

;r  tO’)  •  H,,  ,  J  ( r )  =  e  1,1  a  •  o  »1*'1 

C  ‘  1  <  t  A  ♦  -l 


1,k,l 


K 


L 


i' 

1  io 


(f)  “  n  c 


-  i'  i  A  t 


k,  1 


k-1 


3  =  1 


L 


w-L  E 


k-1 


1=1 


r.vkil(f)  o 


iWiH 


k,  1 


Here  At.  ,  *-  1  o  i,  1-t1  ,  ,  is  the  time  reouired 

r ,  l  < , 1  i.  if  K|  i 

for  the  wavefront  with  the  direction  index  k  end  vo'locit 
index  1  to  neon  f tom  sensor  1  to  sensor  2.  The  expreo- 
s  j on  for  the  coherence  in 


y 


/ri2^f  ^  P?1  ^  \  1/? 

(r”"CfT  v^Xry 


and  accord inr  to  this  for 


inula  and  the  forc/;oin.T  the  coherence  estimate  will  be 
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K  I.  K  L 


DATA  AND  R38ULTS 

In  the  follow  in.-:  we  chall  compare  some  of  the  theo¬ 
retically  computed  coherence  curves  with  coherence  esti¬ 
mates  of  re  1  data  collected  at  Oyer  subarv.uy  (Fir;.  2). 
The  theoretical,  continuous  noi.ee  field  has  been  approxi¬ 
mated  by  3b  discrete  azimuth  directions:  r  nrt  7  individual 
wave  trains  associated  with  each  direction.  The  nowc’ 
densities  of  the  wave trains  decrease  exponentially  from 
a  mximuia  value  at  a  velocity  of  3.3  km/s.  The  total  ve¬ 
locity  span  is  3. 0-3. 6  km/s.  (Fir.  1).  The  cross  and 
auto  spectral  power  estimates  oi  the  r<  J.  data  have  boon 
obtained  by  averr.^in/:  over  50  nonoverlsppinr;  blocks, 
each  of  51.2  see.  Thus  o  ch  estimate  covers  4?  2/3  min 
of  recording  startinr  at  the  times  pivon  on  the  firurc'-. 
The  sampling  interval  was  0.05  see. 


1/2 
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Fi<r  1  shows  a  weather  map  for  May  18,  1968. 
Meteorologically  this  is  a  very  quiet  day  and  by  cx- 
pcriencc  we  do  not  expect  a  very  anisotropic  noise 
field  duo  to  Atlantic  or  coastal  sources  under  such 
conditions.  An  assumed  rower  distribution  at  Oyer 
is  shown  on  too  of  Pi;:.  4.  Even  if  we  assume  that  the 
noise  field  has  a  maximum  in  one  direction,  there  is 
no  reason  to  believe  that  there  is  a  minimum  in  the 
opposite  direction,  'f'herefore  wc  propose  isotropic 
condition  around  the  opposite  half  periphery.  As  wj 
ace  there  is  <a  very  sjood  fit  between  the  experiment.?! 
data  for  1 A 1  -  10Y  and  the  theoretical  curve  calculated 
for  a  sensor  combination  pointing  towards  the  maximum 
noise  power.  The  result  is  a  support  for  the  noise  model 
proposed  and  if  via  use  this  model  it  should  be  located 
with  its  maximum  in  a  north-rest  direction. 

In  the  follo./in."  we  present  examples  of  a  more  ani- 
soti  pic  model  and  for  conra: -iron  vie  have  selected  co¬ 
herence  estimates  made  on  a  da  •  wiih  dominatin':  Atlantic 
and  coastal  noise  sources  (f.i - .  y).  In  the  upper  part 
of  Fiy;.  6  is  /riven  a  proposed  newer  distribution  for  the 
weather  situation  shown  in  hi;;.  5.  The  coherence  curves 


of  Firs.  0  and  7  demonstrate  the  general  incr  asc  in 


coherence  connarcd  to  the  more  isotronie  situation. 

This  is  particularly  pronounced  when  the  sensor  pair 
is  oriented  towards  the  maximum  noise  power,  jf  we 
look  at  a  sensor  combination  abreast,  of  the  maximum 
noise  concentration  (Fie;.  7),  the  coherence  drops  ra¬ 
pidly  from  a  hiyh  start  value. 

One  may  object  that  the  fit  between  the  real  esti¬ 
mates  and  the  theoretical  curves  in  Firs.  6  ana  7  is 
not  verv  rood.  jt  should  bo  nointed  out  that  there 
has  been  no  attempt  made  to  ret  a  better  fit,  for  in¬ 
stance  by  vary in r  details  in  the  parameters.  This  has 
several  reasons:  The  noise  fields  an.  estimated  the 
indirect  way  —  through  the  coherence  —  and  for  comou— 
tat  ion'  1  reasons  the  power  fields  were  defined  usiny 
assumptions  which  are  not  valid  in  general  (c.r.  white¬ 
ness).  Furthermore,  finer  details  in  the  noise  field 
structure  can  not  bo  explored  when  one  is  usinr  only 
two  or  three  space  points. 
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In  this  paper  we  have  used  a  most  co  arnonly  mea~ 

SUrea  °arameter  “  tho  coherence  -  to  estimate  the  noise 
power  distribution.  The  procedure  applied  only  allowed 
tho  corroboration  of  noire  models,  which  were  crude 
approximations  to  the  actual  noise  fields. 

In  the  way  suggested  the  noise  field  can  be  roughly 
estimated  even  if  the  number  of  space  samnlina  points 

-mall  (,.  3)  •  nd  one  gets  a  direct  measure  of  the 
noise  anisotropy. 
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9:  "Spectral 
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tf  1 

•  I  • 


The  theoretic:!  Koine-  Hod el. 


Fir.  4. 


Corr,T)r;rir.oa  between  theoretical  coherence  curves 
(oO lie  lira's )  an--  a  real  coherence  estimate  (do*-* 
cc  line).  Vhc  theoretical  curves  have  been  cal¬ 
culated  us  in."  the  azimuthal  power  distribution 
oh  own  on  ton.  0  refers  to  the  curve  calculated 
lor  a  sensor  combination  oointim  towards  the 
maximum  noise  power,  while  90  refers  to  a  senso- 

L°S510?v8t  90  to  this  direction.  In 

boohcioesthe  sensor  separation  was  3  km. 


4  IS 


retical  curve  (solid  lino 
gives  the  coherence  between  two  censors  2. 
apart  and  whose  connection  line  is  inclined  25  de¬ 
grees  relative  to  the  maximum  noise  power  direction 
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Fig.  7.  The  coherence  estimate  between  IFi  and  1F4  compared 
with  a  theoretical  curve  for  two  censors  2,3  km 
npart  and.  making  an  angle  of  80  degrees  with  the 
direction  of  maximum  noise  cower.  The  power  distri¬ 
bution  is  the  same  as  in  Fig.  6. 
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